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Abstract 

We discuss the momentum distributions of inclusively produced particles 
in e^e~ annihilation. We show that the dependence of the position of the 
maxima of the = ln(l/z) spectra on the mass of the produced particles 
follows naturally from the general definition of fragmentation functions when 
energy-momentum conservation is correctly incorporated. 
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e~^e~ annihilation provides an excellent opportunity to study the fragmentation of quarks 
into hadrons. In particular, inclusive measurements of particle spectra allow us to extract 
fragmentation functions from such experiments and to test different theoretical models of 
fragmentation. Until now, fragmentation functions have not been calculated from first prin- 
ciples rather they have to be modeled in some way. Most of the current approaches use 
different algorithms, such as string and shower algorithms, and model fragmentation of a 
high energy quark in two phases, one of which is purely perturbative, describing the radiation 
and branching of the initial quarks and the other describing the subsequent non-perturbative 
hadronisation of the low energy quarks. Here, we follow a different approach. 

Starting from the general definition of the fragmentation functions, and explicitly guar- 
anteeing energy-momentum conservation we discuss the following interesting property of 
inclusive particle spectra in e^e^ annihilation. When the number of inclusively produced 
particles is plotted as a function of ^ = ln{l/z), (where z is the momentum fraction of 
the fragmenting quark carried by the produced hadron) it exhibits an approximate Gaus- 
sian shape around a maximum, The position of the maximum depends both on the 
total centre of mass energy and on the mass of the produced particle While the 

shape and the energy dependence of the spectrum can be understood in perturbative QCD, 
as a consequence of the coherence of gluon radiation [^], the position of the maximum is 
a free parameter which has to be extracted from experiment. Our main purpose in this 
paper is to show that the dependence of location of the maximum on the mass of the pro- 
duced particle follows naturally from the general definition of fragmentation functions when 
energy-momentum conservation is correctly incorporated. 

Our starting point is the general definition of fragmentation functions 

l^.iz) = -JZj ^e-'^^«-/^Tr{7+(0|^(0)|Pn;pp„)(Pn;pp„)|^(r)|0)}. (1) 

(Here, we discuss only the twist two part of the unpolarized fragmentation functions.) 7+ 
is defined as 7^^ = -^{l^ + 7^) and the plus components of the momenta are defined as 
p+ = -^{jp + p^). p and Pn are the momenta of the produced particle, P, and associated 
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hadronic system n. Using translational invariance to remove the ^~ dependence in the second 
the matrix element and the integral representation of the delta function and projecting out 
the plus components of ip, we obtain 

-D,{z) = ^ E - - P^] I (0|^+(0) |Pn; pp^)\'. (2) 

Here, the plus projection is defined as ip+ = |7+7_'?/'. Using Eq. (|^) rather than Eq. ([^) has 
the advantage that energy-momentum conservation is built in before any approximation is 
made for the states in the matrix element. This is similar to the case of quark distributions 
where the corresponding expression ensures correct support of the distribution functions as 
discussed in Ref. 0. 

While the detailed structure of the fragmentation function depends on the exact form of 
the matrix element some general properties follow already from Eq. (|^). The delta function, 
for example, implies that the function, Dq{z)/z, peaks at 

Here, M and M„ are the mass of the produced particle and the produced system, n, and we 
work in the rest frame of the produced particle. Here, we consider the intermediate state as a 
state having a definite mass. (In general, we have to integrate over a spectrum of all possible 
masses.) We see that the location of the maxima of the fragmentation function depends on 
the mass of the system n. While the high z region is dominated by the fragmentation of a 
quark into the final particle and a small mass system, large mass systems contribute to the 
fragmentation at lower z values. 

We can go a step further and eliminate the (5-function in Eq. by integrating over the 
momentum of the unobserved state n. 

1 1 

-Dq{z) = — = / MPnK0|^+(0)|Pn;pp„)|2, (4) 

Z Zy Z J Pmin 

with 

_ M\l-zf-z-Ml 

' 2Mz{l-z) ^> 
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The significance of Eqs. (^) and (|5D is that Dq{z) vanishes for both z ^ 1 and z ^ 0. Thus, 
fragmentation functions have the correct support. It is interesting to see how the integration 
region depends on the momentum fraction, z, for various values of the produced particle 
system, M„. In Fig. |1], we plot Pmin for the production of protons as a function of z for 
different values of M„. Pmin = gives the value of z at which the contribution of a given 
Mn to the fragmentation function is largest. This gives Zmax according to Eq. Further, 
the region of z where the lower integration limit is sufficiently small that the integral will 
be significant, becomes narrower with increasing M„. Thus, large mass states contribute to 
the fragmentation function at low z and only in a very narrow region of z. 

At a given centre of mass energy, there will be a maximum value for the mass of the 
intermediate state which can be produced in the fragmentation. This maximal mass deter- 
mines the "lower" edge of the spectrum. Q We can use Eq. (j^) to estimate the maxima 
of the ^-distribution associated with this particular mass. This maximum determines the 
maximum of the fragmentation function in first approximation. Although M„ is not known, 
it should be proportional to the available total energy Ecu- However, the precise value of 
Mn is not needed if we are only interested in the relative position of the maxima of the ^ 
distribution of different particles. Taking the difference of the maxima of the ^ = ln(l/^;) 

^ Since the ^-distribution is given by da/d^ = zdajdz ~ zD(z) it is proportional to times 
D{z)/z. Although Eq.(^) describes the location of the maximum of the distribution D{z)/z we can 
expect that Eq.(^ is also a good approximation for the ^-distribution, since the z region where the 
lower integration limit (pmin), is sufficient small, is very narrow for large masses, M„. Thus, the 
contributions from a given M„ to the fragmentation functions are very narrow functions in z for 
large M„. Note, that the square of the matrix element in Eq. (^) must decrease faster then l/p^ 
in order to give finite result for the fragmentation functions. Eq. (^) gives the maximum of the 
^-distribution exactly in the limiting case when the contribution of a given Mn to the fragmentation 
function is a 5-function. 
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distributions of two different particles, a and b, the dependence on the unknown value of 
Mn drops out for sufficiently large M„. It follows from Eq.(|]) that 

Ar ~ In (Ml±^] + In ^ ~ In ^ (6) 

Thus, the difference of the maxima is determined by the logarithm of the ratio of the masses 
of the produced particles. Since the value of M„ for finite centre of mass energies are in 
general different for mesons and baryons, Eq. (P) will be only valid for the difference of the 
maxima of the mesons or baryons separately. We calculated the maxima of the ^ distributions 
using this formula and using the maxima of the 1]' and that of the proton distributions as 
a reference value for mesons and baryons, respectively. The results are compared to the 
experimental data in Fig. |^. The location of the maxima as a function of the mass of 
the produced particles is reasonably well described both for mesons and baryons. 

In conclusion, we have shown that the dependence of the position of the maxima of the 
^ = ln(l/;z) spectra on the mass of the produced particles follows from the general definition 
of the fragmentation functions and from energy-momentum conservation. Our results are in 
remarkably good agreement with the data. 
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FIG. 2. Location of the maxima of the distributions as a function of the particle mass. The 
full and open symbols represent mesons and baryons, respectively. The data are from Refs. 
The solid and dashed lines are the predictions of Eq. for mesons and baryons adjusting the 
normalization to and to .^*, respectively. 
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